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We measure the mass of the top quark using top quark pair candidate events in the lepton+jets 
channel from data corresponding to 1 fb _1 of integrated luminosity collected by the DO experiment 
at the Fermilab Tevatron collider. We use a likelihood technique that reduces the jet energy scale 
uncertainty by combining an in-situ jet energy calibration with the independent constraint on the 
jet energy scale (JES) from the calibration derived using photon+jets and dijet samples. We find 
the mass of the top quark to be 171.5 ± 1.8(stat.+JES) ± l.l(syst.) GeV. 
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Since the discovery of the top quark in 1995 [l], a sub- 
stantial effort has gone into measuring and understanding 
its properties. Its large mass suggests a unique role in the 
mechanism of electroweak symmetry breaking. Through 
radiative corrections, a precise measurement of the top 
quark mass, together with that of the W boson, allows 
indirect constraints to be placed on the mass of the stan- 
dard model Higgs boson [2| . A precise knowledge of the 
top quark mass could also provide a useful constraint to 
possible extensions of the standard model. It is therefore 
of great importance to continue improving measurements 
of the top quark mass [3|, |J] . 

In this Letter, we present the most precise single mea- 
surement of the top quark mass from Run II of the Fer- 
milab Tevatron collider. It uses a matrix element (ME) 
method with an in-situ jet energy calibration based on 
a global factor used to scale all jet energies and thereby 
the invariant mass of the hadronic W boson [3, 5]. This 
mass is constrained to the well known value of 80.4 GeV 
through the Breit-Wigner function for the hadronic W 
boson in the ME for ti production. The jet energy scale 
is further constrained to the standard scale derived from 
photon+jets and dijet samples within its uncertainties 
through the use of a prior probability distribution. This 
analysis is based on data collected by the DO detector [6[ 
from April 2002 to February 2006 comprising 1 fb^ 1 of in- 
tegrated luminosity from pp collisions at yfs = 1.96 TeV. 

The top quark is assumed to always decay into a W bo- 
son and a b quark producing a W + W~bb final state from 
it production. This analysis is based on the lepton+jcts 
channel with one W boson decaying via W — ► Iv and the 
other via W — ► qq'. This channel is characterized by a 
lepton with large transverse momentum (pr), large mo- 
mentum imbalance due to the undetected neutrino (j> T ), 
and four high-px jets. Events are selected for this analy- 
sis by requiring exactly one isolated electron (muon) with 
p T > 20 GeV and \r)\ < 1.1 (\r)\ < 2), j> T > 20 GeV, and 
exactly four jets with pr > 20 GeV and |?7| < 2.5, where 
the pseudorapidity rj = — In [tan(6*/2)], and 9 is the polar 
angle with respect to the proton beam direction. Multijct 
background, typically originating from lepton or jet en- 
ergy mismeasurements, is further suppressed by requiring 
the lepton direction and j) T vector to be separated in az- 
imuth. Jets are reconstructed using a cone algorithm Q 
with radius R = \J {Ay) 2 + (Acj)) 2 — 0.5 where the y 
is the rapidity. Jet energies are corrected to the particle 
level using corrections derived from photon+jet and dijet 
samples. Jets containing a muon are assumed to origi- 
nate from semileptonic b quark decays and corrected by 
the muon momentum and average neutrino energy. At 
least one jet is required to be tagged by a neural-network 
based algorithm [8| as a &-jet candidate. The tagging effi- 
ciency for b jets is ~ 50% with a misidentification rate of 
~ 1% from light jets. A total of 220 events, split equally 
between the electron and muon channels, satisfying these 
criteria is selected. 



The top quark mass is determined from the data sam- 
ple with a likelihood method based on per-event proba- 
bility densities (p.d.'s) constructed from the MEs of the 
processes contributing to the observed events. Assuming 
only two processes, it and W+jets production, the p.d. 
to observe an event with measured variables x is 

P cvt = A(x) [/P sig (x; m t ,k jcs ) + (1 - /)Pbk g (a;; fcjes)] , 

where the top quark mass mt, jet energy scale factor 
kj es dividing the energies of all jets, and observed signal 
fraction / are the parameters to determine from data. 
P s i g and Pbkg are, respectively, p.d.'s for it and FF+jets 
production. Multijet events satisfy Pbkg 3> Pig and are 
also represented by Pbkg- A(x) is a function only of x and 
accounts for the geometrical acceptance and efficiencies. 

P s ig and Pbkg are calculated by integrating over all pos- 
sible parton states leading to the measured set x. In addi- 
tion to the partonic final state described by the variables 
y 7 these states include the initial state partons carrying 
momenta qi and q2 in the colliding p and p. The inte- 
gration involves a convolution of the partonic differential 
cross section da(y; mt) with the p.d.'s for the initial state 
partons f{qi) and the transfer function W(y, x; fcj es ): 



Psig = J ^2da(y 1 m t )dqidq 2 f(qi)f(q2)W(y,x;kj, 



where the sum runs over all possible initial state par- 
ton flavor combinations. f(qi) includes parton den- 
sity functions (PDFs) for finding a parton of a given 
flavor and longitudinal momentum fraction in the p 
or p (CTEQ6L1 and parameterizations of the 

p.d.'s for the transverse components of qi derived from 
PYTHIA [l(| • Jet fragmentation effects and experimental 
resolution are taken into account by W{y, x; kj es ), repre- 
senting the p.d. for the measured set x to have arisen 
from the partonic set y. The normalization factor TV, de- 
fined as the expected observed cross section for a given 
(mt, fcjes); ensures A(x)P S i g (and ultimately P ev t) is nor- 
malized to unity. 

The differential cross section term in P s j g is calculated 
using the leading order ME for qq — > ti. After all en- 
ergy and momentum constraints are taken into account, 
this term is integrated over the energy associated with 
one of the quarks from the hadronic W boson decay, 
the masses of the two W bosons and two top quarks, 
and the energy (l/px) of the electron (muon). It is 
summed over 24 possible jet-parton assignments each 
carrying a 6-jet tagging weight llj and over the neu- 
trino solutions. W(y, x; kj cs ) is the product of five terms 
for the four jets and one charged lepton. The jet terms 
are parameterized in terms of jet energy with a func- 
tion involving the sum of two Gaussians whose param- 
eters depend linearly on parton energy. The term for 
the charged lepton is parameterized as a Gaussian dis- 
tribution in energy (l/px) for electrons (muons). All 
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FIG. 1: Comparison between data and MC 2-jet (m2j) and 
3-jet (m3j) invariant mass distributions. 



parameters for W(y, x; fcj OS ) are derived using fully simu- 
lated Monte Carlo (MC) events. The normalization cross 
section er" bs = j A(x)P s - lg dx = a u (m t ) (A(m t , fcj es )) is 
calculated using the total cross section corresponding to 
the ME used and the mean acceptance for events whose 
dependencies on mt and fcj es are determined from MC 
events. 

The differential cross section term in Pbkg is calculated 
using the W+4 jets MEs provided by VECBOS [IH ■ The 
integration is performed over the energies of the four par- 
tons producing the jets, the W boson mass, and the en- 
ergy (1/pt) of the electron (muon) summing over 24 pos- 
sible jet-parton assignments and two neutrino solutions. 
The transverse momenta of the initial state partons are 
assumed to be zero. 

P s i g and Pbkg are calculated using MC techniques 
on a grid in (mt,fcj es ) having spacings of 1.5 GeV and 
0.015, respectively. At each grid point, a likelihood func- 
tion, L{x\ mt, fcjes, /), is constructed from the product of 
the individual event p.d.'s (P ev t) and / is determined 
by minimizing — lnL at that point. m*, fcj cs ) is 

then projected onto the m t and fcj es axes according to 
L(x;m t ) = / L(x;m t ,k ics )G(k ies )dkj cs and L(x;k jcs ) = 
J L(x; m tl k- ]cs )dm t . The prior G(fcj OS ) is a Gaussian func- 
tion centered at fcj es = 1 with width a = 0.019 deter- 
mined from the photon+jets and dijet samples used in 
the standard jet energy scale calibration. Best estimates 
of mt and fcj os and their uncertainties are then extracted 
from the mean and rms of L(x;mt) and L(x; fcjes). 

The measurement technique described above is cali- 
brated using MC events produced with the ALPGEN event 
generator [13[ employing pythia for parton showering 
and hadronization and implementing the MLM matching 
scheme Il4j ■ All generated events are processed by a full 
GEANT [ToJ detector simulation followed by the same re- 
construction and analysis programs used on data. Fig. [T] 
shows comparisons of the 2-jet and 3-jet invariant mass 
distributions between data and MC using it events with 
a true top quark mass (m}™ ) of 170 GeV. These are 
calculated using jets assigned as the decay products of 
the top quark and W boson from the hadronic branch 




m true_-| 7 0GeV k|es 6 - 1 

FIG. 2: Mean values of mt and fcjes from ensemble tests ver- 
sus true values parameterized by straight lines. Dashed lines 
represent identical fitted and true values. 



in the permutation with the largest weight (defined as 
the product of P s i g and the 6-jet tagging weight) around 
the peak of L(x; mt, fcjes)- MC distributions are normal- 
ized to data distributions with / = 0.74 determined from 
data. The background includes simulated VF+jets events 
and data events selected from a multijet enriched sam- 
ple. The latter comprises 12% of the total background 
based on estimates from data. The estimated number of 
ti events (e + jets: 91 ± 9, /x + jets: 71 ± 8) agrees with 
the expectation (e + jets: 89 ± 6, fi + jets: 73 ± 5). 

Five ti MC samples are generated with m* ruc = 160, 
165, 170, 175, and 180 GeV, with six more produced from 
three of these by scaling all jet energies by ±5%. P s i g and 
Pbkg are calculated for these events from which pseudo- 
experiments fixed to the number of data events are ran- 
domly drawn with a signal fraction fluctuated according 
to a binomial distribution around that determined from 
data. The mean values of m t and fcj OS for 1000 pseudo- 
experiments are shown in Fig. [5] as functions of the true 
values and fitted to a straight line. The average widths of 
the mt and fcj CS pull distributions are 1.0 and 1.1, respec- 
tively. The pull is defined as the deviation of a measure- 
ment from the mean of all measurements divided by the 
uncertainty of the measurement per pseudo-experiment. 
The measured uncertainties in data are corrected by the 
deviation of the average pull width from 1.0. 

L(x;mt) and L{x;k- }es ) for the selected data samples 
are calibrated according to the parameterizations shown 
in Fig. [21 L(x]m t ) is shown in Fig. [31(a) with a mea- 
sured m t = 171.5 ± 1.8(stat.+JES) GeV. The measured 
fcjes = 1.030 ± 0.017 represents a 1.2 a shift from fcj CS = 1 
where a is the sum in quadrature of the width of G(fcj C s) 
and the uncertainty of the measured fcj cs . Fig. [31(b) com- 
pares the measured uncertainty for mt with the expected 
uncertainty distribution from pseudo-experiments in MC 
assuming TO( 1UC = 170 GeV. 

To verify the in-situ jet energy calibration, we repeat 
the analysis on data by fixing fcj os to the measured value 
and removing the W boson mass constraint, replacing 
L(x; m t , fcjes, /) with L(x; m t ,m w , /)• Psig and P b k g are 
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FIG. 3: (a) Projection of data likelihood onto the mt axis with 
best estimate shown, (b) Expected uncertainty distribution 
for mt with measured uncertainty indicated by the arrow. 



TABLE I: Summary of systematic uncertainties (symmetrized 
based on the larger of the two values in each direction) . 



Source 


Uncertainty (GeV) 


Physics modeling: 




Signal modeling 


±0.40 


PDF uncertainty 


±0.14 


Background modeling 


±0.10 


b- fragmentation 


±0.03 


Detector Modeling: 




6/light response ratio 


±0.83 


Jet identification and resolution 


±0.26 


Trigger 


±0.19 


Residual jet energy scale 


±0.10 


Muon resolution 


±0.10 


Method: 




MC calibration 


±0.26 


6-tagging efficiency 


±0.15 


Multijet contamination 


±0.14 


Signal contamination 


±0.13 


Signal fraction 


±0.09 


Total 


±1.07 



now calculated on a grid in (mt,mw) having spacings of 
1.5 GeV and 1 GeV, respectively. L(x; m t ) and L(x; mw) 
are calculated in the same way as for the grid in (mt,fcj es ) 
except that no prior probability distribution is used for 
L(x;mt). We find mw = 80.3 ± 1.0 GeV which is consis- 
tent with the constraint of 80.4 GeV. 

Systematic uncertainties are evaluated for three cate- 
gories. The first category involves the modeling of MC 
tt and VF±jets events and includes uncertainties in the 
modeling of extra jets due to radiation in tt events, 
the distribution shapes and the heavy flavor fraction in 
VF±jets events, b fragmentation, and the PDFs used in 
generating events. The second category is associated 
with the simulation of detector response and includes 
possible effects due to the energy and \n\ dependence 
of the jet energy scale unaccounted for by the in-situ 
calibration, uncertainties in the modeling of the relative 
calorimeter response to b and light quark jets, and un- 
certainties associated with the simulation of jet energy 



resolution and reconstruction efficiency, muon px resolu- 
tion, and trigger efficiency. The third category is related 
to assumptions made in the method and uncertainties in 
the calibration and includes possible effects due to the ex- 
clusion of multijet events and non-lepton±jets tt events 
from the calibration procedure, uncertainties in the signal 
fraction used in ensemble tests, and uncertainties associ- 
ated with the parameters defining the calibration curve. 
Contributions from all these sources are summarized in 
Table Hand sum in quadrature to ±1.1 GeV. 

The leading sources of uncertainty in Table U are those 
associated with the 6/light response ratio and signal mod- 
eling. The first of these is evaluated by estimating the 
possible difference in this ratio between data and MC 
and scaling the energies of all jets matched to b quarks 
in a MC tt sample by this amount. The analysis is re- 
peated for this sample and the difference in nit from that 
of the unsealed sample taken as the uncertainty. The 
uncertainty associated with the modeling of additional 
jets in tt events is evaluated using both data and MC 
samples. Using MC tt events, the fraction of tt signal 
events with > 5 jets is varied such that the ratio of 4-jet 
to > 5-jet events in MC matches that in data including 
its uncertainties. The difference in the resulting m t from 
that of the default sample is taken as the uncertainty. 
Using data, this is done through ensemble tests in which 
a fixed number of > 5-jet events not used in the mea- 
surement are randomly drawn for each experiment and 
combined with the default sample of 4-jet events. The 
ensemble tests are repeated for different fractions of > 5- 
jet events constituting up to 30% of each experiment. m t 
for the default sample is compared with the mean from 
each ensemble test and the largest difference taken as the 
systematic uncertainty. Both procedures yield consistent 
estimates for this systematic uncertainty. 

In summary, we present a measurement of the top 
quark mass using it lepton±jets events from 1 fb -1 of 
data collected by the DO experiment. Using a ME tech- 
nique combining an in-situ calibration of the jet energy 
scale with the calibration based on the photon±jets and 
dijet samples gives us 

m t = 171.5 ± 1.8(stat.±JES) ± l.l(syst.) GeV, 

representing the most precise single measurement to date. 
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